iT邦幫忙

第 11 屆 iThome 鐵人賽

DAY 16
0
Software Development

開源的GIS實作系列 第 16

[day-16] 影像處理 - 空間濾波器

  • 分享至 

  • xImage
  •  

前言

我們接續執行完ODM產出成果的進度,將資料爬出來來並且處理。。

檔案結構

main.py
README.md
upload_file/
results/
├── odm_dem/  # DSM放置處
└── odm_orthophoto/ #正射鑲嵌影像
results_with_dtm/

影像的空間濾波

影像需求

空間濾波器可以強調影像中的某些特徵或模糊影像,很多的套件中都也相關的資料濾波器可以使用,但是有一點要注意的是,濾波完的影像要寫回去原來的檔案,並且需要存有metadata這是比較需要注意的地方。

variance filter

Variance 的公式

理論公式如下,可以凸顯影像的邊界。
https://chart.googleapis.com/chart?cht=tx&chl=Variance%20%3D%20E(X%5E2)%20-%20%20E(X)%5E2

實作code

def std_convoluted(image, N):
    im = np.array(image, dtype=float)
    im2 = im**2
    ones = np.ones(im.shape)

    kernel = np.ones((2*N+1, 2*N+1))
    s = scipy.signal.convolve2d(im, kernel, mode="same")
    s2 = scipy.signal.convolve2d(im2, kernel, mode="same")
    ns = scipy.signal.convolve2d(ones, kernel, mode="same")

    return np.sqrt((s2 - s**2 / ns) / ns)

程式碼來自:https://stackoverflow.com/questions/25910050/perform-local-standard-deviation-in-python


上一篇
[day-15] import 自己的函式庫
下一篇
[day-17] 影像處理串接
系列文
開源的GIS實作30
圖片
  直播研討會
圖片
{{ item.channelVendor }} {{ item.webinarstarted }} |
{{ formatDate(item.duration) }}
直播中

尚未有邦友留言

立即登入留言